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We discuss anarchy in the type-I seesaw framework with two right-handed neutrinos. For a 
seesaw scale around eV, the new fermions make up two light sterile neutrinos to accommodate 
various experimental anomalies. The distributions of mixing angles and masses according to anarchy 
with naturai measures are in agreement with global fits. The framework prediets the absence of 
neutrinoless doublé beta decay and one vanishing neutrino mass, and can therefore be tested in 
future experiments. 
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I. INTRODUCTION 



Neutrino oscillations have established the need for 
massive neutrinos and hence physics beyond the Stan- 
dard Model (SM). While most of the data can bc 
explained in a minimal framework of three active 
neutrinos — realized for example with a type-I see- 
saw mechanism [l[ — some experiments and observations 
strongly favor one or two additional light neutrinos, with 
masses around eV and 0(0.1) mixing Those hints 

stoni from completcly independent ficlds, namely parti- 
cle physics, cosmology and astrophysics. We refer to the 
exhaustive overview of the situation provided by Ref . Q , 
a White Paper dedicated to the subject. Since these new 
states are not allowed to contribute to the Z-width, they 
have been dubbed sterile neutrinos. A large number of 
experiments will be coming up in the months and years to 
come, with the aim to investigate their possible presence 
with a variety of methods and approaches [f| . Obviously, 
proving the existence of light sterile neutrinos would bc 
a sensational and groundbreaking discovery. 

Sterile neutrinos aside, the mixing of active neutrinos 
has been of much interest to the physics community, as 
the peculiar values for the mixing angles seem to hint at 
an underlying symmetry principle. The astonishing pre- 
cision of neutrino experiments nowadays has however put 
a slight dent in this idea, as models seem to get increas- 
ingly complicated in order to be valid. On the other end 
of the model-building spectrum, the complete absence 
of a symmetry behind lepton mixing has also been pro- 
posed both for an effective theory of neutrino mass 

and for the type-I seesaw case. This "anarchy solution" 
(henceforth called active anarchy) has been successful, 
especially in light of the rather large reactor mixing an- 
gle 6<i3 [§|. 

Taking the various hints for sterile neutrinos seriously, 
it is of obvious interest to try to accommodate these 
states in models. While fiavor symmetry models that 
include sterile neutrinos exist Q, the counterframework 



of anarchy has not yet been discussed in the context of 
light sterile neutrinos in detail. This will be amended by 
the letter at hand. 

Instead of simply adding two light sterile neutrinos to 
the three active ones, we will work in a much more eco- 
nomical scenario, namely by assuming that these states 
are the right-handed neutrinos of the type-I seesaw mech- 
anism (eV-seesaw 0]). The most minimal case which can 
accommodate the active neutrino data is then when only 
two such states are added to the picture.^ This implics 
two light sterile neutrinos, one masslcss active neutrino, 
no neutrinoless doublé beta decay and only a small effec- 
tive electron neutrino mass. A phenomenological study of 
this straightforward scenario has recently been provided 
in Ref. [11| ■ We will study the implications of anarchy in 
this framework and present the statistical distributions 
for ali observables. We show that the large active-sterile 
mixing angles needed for various anomalies can be natu- 
rally obtained. Thus, the attractive scenario of anarchy 
can also be extended to the case of light sterile neutrinos. 



II. ANARCHY 

Two right-handed singlet neutrinos Ai. 2 modify the 
SM-Lagrangian by the following terms 
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where sums over j, k = 1,2 and a = e, fi, r are under- 
stood. After electroweak symmetry breaking, using the 
usuai SM Higgs doublet H — > (0, v) T ~ (0,174GcV) T , 
we arri ve at the 5x5 Majorana mass matrix for the neu- 
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trai fermions in the basis {vi : N c ) T 

Mf « n - [vVnDyVl Mn ) (2) 
= U* diag(mi, rri2, rri3, m^, m.5) W . 

Here, we used the singular value decomposition for the 
Yukawa coupling matrix 

^ = ^(o l t)Vl = VnDrVL (3) 

and introduced a unitary matrix U to diagonalize A4f u ii. 

In the anarchy framework we assume that the unitary 
3 x 3 (2 x 2) matrix Vi (Vr) is distributed according 
to the Haar measure of U(2>) (U(2)), which can then be 
compared to the experimental values. § As far as the sin- 
gular values yj and eigenvalues of Mr go, we assume a 
distribution according to the linear measures as derived 
in Refs. 0, EU 

dM R oc \M? - M$\MiM 2 dAfidAf 2 , 
dD Y oc {y{ - ylfylvl à yi dy 2 , 

with the boundary conditions tifòyDy) — J2ìVì — 
y 2 and tr (M R M R ) = J2 t M ì < M o- A survey of 
more complicated measures including weighting func- 
tions, analogous to Ref. [l2[ for active anarchy, lies bc- 
yond the scope of this letter. 

Note that this is the proper way of defining anarchy 
in the case of two right-handed neutrinos, even outsidc 
of the seesaw limit M.r S> vY v . Taking for example 
the full mass matrix A^f u ii to be distributed according to 
1/(5) (as in Ref. [l3j]) would make Majorana masses for 
the vi necessary, for which we would have to introduce 
a scalar triplct (type-II seesaw), complicating the study 
considerably. Additional mechanisms that decouple some 
heavy right-handed neutrinos or impose anarchy only in 
a subsector could also be constructed [l4|, but we will 
only consider the minimal framework here. 

In Eq. (T2J) we have mi = 0, and the only free param- 
eters are the scales yo and M . M can be fixed to give 
niQ ~ 1 eV, while the Yukawa scale yo is fixed to render 
777,3 — 0.05 cV. The other two masses and ali the mixing 
anglcs and phases are then distributed in a known way 
and can be compared to measurements. Defining the pa- 
ramctcr e = vyo/Mo, we can estimate the expected values 
of the observables 

m 2 ,3 - Mo e 2 (1 + 0{e 2 )) , 7774,5 ~ M (1 + 0{e 2 )) , 
U ai , U a5 ~ e (1 + 0(e 2 )) , [/pmns - (1 + 0(e 2 )) . 
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Figure 1: Distribution of mass ratios Rij = Am^/ Am^ for 
the values M = 1.6 eV, y = 1.5 x 10~ 12 . The red/dotted 
(red/dot-dashed) distribution corresponds to R23 in the see- 
saw limit with two (three) sterile neutrinos. Vertical lines 
denote best-fìt values, shaded areas 90% C.L. (sterile) and 
3<r (active) ranges. 
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Figure 2: Distribution of masses for the same values as in 

Fig.ru 

Here, J/pmns denotes the upper left 3x3 submatrix 
of U (corresponding to the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) mixing matrix), which contains the ac- 
tive mixing angles. In the eV-seesaw limit that we are 
intcrcstcd in. e < 0.2. The caso 5 — » is the canonica! 
seesaw limit. We can estimate that in case of eV-seesaw 
the ratios 

i? 23 = A7772 1 /A?77 3 ! 1 , R 45 = Am^/Am^ , (5) 

will only receive percent corrections compared to the see- 
saw limit. This can be seen in Fig. [TJ where the distri- 
butions for R23 are shown for e ~ 0.2, e <C 0.2 (77 = 2) 
and e <§; 0.2 (77 = 3): the hierarchy of i?23("- = 2) is 
pulled closer to R2z{n = 3) due to the next-to-leading 
order seesaw contributions.^ 

As interesting numerical values throughout this letter 
we use M = 1.6 eV, y = 1.5 x 10~ 12 (i.c. e ~ 0.16). 



§ Note that it does not matter whethcr we takc a diagonal Mn = 
diag(Mi,M 2 ) or VMrV t with Haar-distributed V G U(2), be- 
cause our framework is by construction basis-independent. 



■ Here and in the following, ali shown distributions are properly 
normalized. The sampling procedure follows Ref. lidi . 
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In our Figs. HHH we show the distributions of some ob- 
scrvablc quantities as well as currcnt global-fit values. 
The best-fit values and 3(7 ranges for the active neu- 
trino parameters are taken from Ref. [l5j . the values for 
the active-sterile mixing elements | Cq.4,5 | from Ref. [Jj|, 
and the masses Am^, Am^ from Ref. Q (Fig. 71). 
Using other global fits for these parameters, e.g. from 
Refs. Hi], docs of course not change the qualitative dis- 
cussion of this letter. We further note that the precision 
in the parameters associated to the sterile neutrinos is 
much weaker than the one for the active ones. 

Let us emphasize that while the best-fit value for R45 
is somewhat larger than the anarchy prediction, this ten- 
sion is alleviated once we consider e.g. the 90% C.L. con- 
tour. For example, there seems to be an interesting re- 
gion in parameter space Am^ ~ lcV 2 , Am^ ~ 6cV 2 [H 
that gives roughly R 45 =0.17, perfect for sterile anarchy. 
Obviously, the agreement with standard ACDM cosmol- 
ogy worsens in this case due to the large sum of neutrino 



masses 



17j . Note also that the preference of two over one 
sterile neutrino has weakened |18( due to a recent update 
from the MiniBooNE experiment [l9| , which reduces the 



previous tension between the v„ 



and v„ 



data. 



Furthcrmore, a stringent 99% C.L. limit for the active- 
sterile mixing angle of sin 6* < 0.07 in the relevant mass 
rango was given very recently by the ICARUS experi- 
ment [20|. However, from Fig. [3] (bottom) we see that a 
smaller mixing angle is also fine in sterile anarchy. 

Let us discuss the distribution of the mixing matrix U . 
Since we are not working in the seesaw limit, or in an 
effective theory, the upper left 3x3 submatrix C/pmns is 
not distributed according to the Haar measure of U(3) 
as in active anarchy Q 



dC/pMNS dsT2 d s 23 de* 3 d5 da d/3 . 



(6) 



(which would imply the same distribution for ali |ì7 a j|, 
a = e,fi,T, j = 1,2,3), but shows small deviations. 
From the unitarity of U — namely J^- |£/ Q i| 2 = 1 — we ex- 
pect diminished values for |£/ Q i,2,3| upon increasing the 
active-sterile mixing |{/ a 4,s|' The diagonalization condi- 
tion ^ U^rrii = — related to the upper left 3x3 zero 
matrix in *Mf u u — then shows that mainly \U a 3\ is sup- 
pressed (due to m 3 > m 2 ). This is illustrated in Fig. [3] 
(top), where we see that \Ues\ (and therefore #12) is stili 
approximately distributed likc in the seesaw limit, while 
\U a s\ is drawn to slightly smaller values. This helps the 
agreement of anarchy with the rather small (in compar- 
ison to #12,23) reactor mixing angle \U e z\ — S13 — 0.16, 
but makes the agreement with the (comparably impre- 
cise) atmospheric mixing angle #23 a little worse, both 
however insignificant. 

The active-sterile mixing elements |[/ Q 4.5| are large, as 
expected from the scaling |J7 ct 4 i 5| ~ e ~ 0.1 (bottom 
Fig. [3]). The suppression of \U a s\ compared to \U a 4\ is 
due to the mass hierarchy in Mr, as the entries scale 
with and l/^/mJ, respectively @. We are led to 

conclude that sterile anarchy inherits the success of active 
anarchy for the PMNS mixing matrix, and is further fully 
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Figure 3: Distribution of some mixing elements for the val- 
ues Mo = 1.6 eV, yo = 1.5 x 1CP 12 . The vertical lines rep- 
resent best-fit values, the shaded areas are 3<r ranges (top 
from Ref. [TU, bottom for a — fi from Ref. [JJ|). In the see- 
saw limit, U a 5, U a 4, go to zero, while U a 3 and U a 2 converge 
roughly to the distribution for U e 2 (top). 



compatiblc with large active-sterile mixing anglcs. 

Having discussed the distributions of masses and mix- 
ing angles, we now bricfly turn to other observables. First 
of ali, as in any anarchy scenario, the normal mass or- 
dering is preferred. In our case only about 5 % of the 
cases give the inverted ordering. The relevant phase 
for short-baseline experiments arg(C/* 4 C/e5c7 M 4t/* 5 ) is dis- 
tributed uniformly from zero to 27r — as expected from 
phases in anarchy — and can therefore not be used to test 
this framework. The rate of neutrinoless doublé beta 
decay [UJ vanishes; this is because ali neutrino masses 
are far below the momentum exchange q 2 ~ (100 MeV) 2 
relevant for the decay, so the matrix element will be pro- 
portional to (M.f u u)ee = 0. Kurie-plot (beta decay) ex- 
periments can test the effective electron neutrino mass 



| 2 m 2 



(7) 



with a current upper limit of 2.3 eV at 95% C.L. [22| . 
The distribution of this quantity is shown in Fig. [H so 
mp is expected to be rather small; correspondingly, a 
discovery in the KATRIN experiment [23| (detection po- 
tcntial m/3 = 0.35 eV) would pretty much exclude sterile 
anarchy. 
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Figure 4: Distribution of mp for the values Mo = 1.6 eV, 
yo = 1.5 x 10~ 12 . The shaded region corresponds to the 
90% C.L. design sensitivity of KATRIN [13|. 



this framework to a minimal scherno with light sterile 
neutrinos. The mass hierarchies and mixing angles are 
surprisingly close to the ones needed to explain a number 
of neutrino anomalies, if the only two parameters of the 
model — the Yukawa and seesaw scales — are chosen ap- 
propriatela Statistical improvements for the sterile neu- 
trino parameters are necessary to properly evaluate the 
validity of the approach discussed here. Whilc currently 
not really predictive due to the low precision of sterile 
neutrino parameters, the model can be easily excludcd, 
as it hinges on the absence of neutrinoless doublé beta 
decay, no discovery at KATRIN and the existence of two 
light sterile neutrinos, ali of which are falsifiable in the 
near future. This is to be compared to active anarchy, 
which is notoriously hard to test for. 



III. CONCLUSION 

Anarchy for active neutrinos with and without a seesaw 
mechanism has been shown to be in good agreement with 
the measured neutrino mixing parameters. In this let- 
ter we have demonstrated that an eV-scale type-I seesaw 
mechanism with two right-handed neutrinos can extend 
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